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The production of charmed and beauty hadrons in proton-proton and proton-antiproton 
collisions at high energies are analyzed within the modified quark-gluon string model 
(QGSM) including the internal motion of quarks in colliding hadrons. We present some 
predictions for the future experiments on the beauty baryon production in pp collisions 
at LHC energies. This analysis allows us to find interesting information on the Regge 
trajectories of the heavy (66) mesons and the sea beauty quark distributions in the proton. 



1 Introduction 

Various approaches of perturbative QCD including the next-to-leading order calculations (NLO 
QCD) have been applied to construct distributions of quarks in a proton. The theoretical 
analysis of the lepton deep inelastic scattering (DIS) off protons and nuclei provides rather 
realistic information on the distribution of light quarks like u, d, s in a proton. However, to find 
a reliable distribution of heavy quarks like c(c) and especially b(b) in a proton describing the 
experimental data on the DIS is a non-trivial task. It is mainly due to small values of D and B 
meson yields in the DIS at existing energies. Even at the Tevatron energies the B- meson yield 
is not so large. At LHC energies the multiplicity of these mesons produced in pp collisions will 
be significantly larger. Therefore one can try to extract a new information on the distribution 
of these heavy quarks in a proton. In this paper we suggest to study the distribution of heavy 
quarks like c(c) and b(b) in a proton from the analysis of the future LHC experimental data. 

The multiple hadron production in hadron-nucleon collisions at high energies and large 
transfers is usually analyzed within the hard parton scattering model (HPSM) suggested in [HE] ■ 
This model was applied to the charmed meson production both in proton-proton and meson- 
proton interactions at high energies, see for example [3]. The HPSM is significantly improved by 
applying the QCD parton approach [U[S], see details in [Sj and references therein. Unfortunately 
the QCD including the next-to-leading order (NLO) has some uncertainties related to the 
renormalization parameters especially at small transverse momenta pt [6] . 

In [6J [7] we studied the charmed and beauty meson production in pp and pp collisions at 
high energies within the QGSM [5] or the dual parton model (DPM) [5] based on the 1/N 
expansion in QCD [TUIITT]. It was shown that this approach can be applied rather successfully 
at not very large values oipt- In this paper we investigate the open charm and beauty baryon 
production in pp collisions at LHC energies and very small pt within the QGSM to find new 
information on the Regge trajectories of the heavy (cc) and (bb) mesons and the sea beauty 
quark distributions in the proton. 
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2 General formalism for hadron production in pp collision 
within QGSM 

Let us present briefly the scheme of the analysis of the hadron production in the pp collisions 
within the QGSM including the transverse motion of quarks and diquarks in colliding protons 
|12j . As is known, the cylinder type graphs for the pp collision presented in Fig. 1 make the main 
contribution to this process [8j. The left diagram of Fig.l, the so-called one-cylinder graph, 
corresponds to the case where two colorless strings are formed between the quark/diquark 
(q/qq) and the diquark/quark (qq/q) in colliding protons; then, after their breakup, qq pairs 
are created and fragmcntated to a hadron, for example, D meson. The right diagram of Fig.l, 
the so-called multicylindcr graph, corresponds to creation of the same two colorless strings and 
many strings between sea quarks/antiquarks q/q and sea antiquarks/quarks q/q in the colliding 
protons. The general form for the invariant inclusive hadron spectrum within the QGSM is 
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Figure 1: The one-cylinder graph (left diagram) and the multicylindcr graph (right diagram) 
for the inclusive pp — > hX process. 
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where E, p are the energy and the three-momentum of the produced hadron h in the laboratory 
system (l.s.) of colliding protons; E* , s are the energy of h and the square of the initial energy in 
the c.m.s of pp; x,pt are the Fcynman variable and the transverse momentum of h; a n is the cross 
section for production of the n-Pomcron chain (or In quark-antiquark strings) decaying into 
hadrons, calculated within the "eikonal approximation" [14] . Actually, the function 4> n (x,pt) 
is the convolution of the quark (diquark) distributions in the proton and their fragmentation 
functions (FF) , see details in [51 O El H2]- To calculate the interaction function <f) n (x,pt) 
we have to know all the quark (diquak) distribution functions in the nth Pomeron chain and 
the FF. They are constructed within the QGSM using the knowledge of the secondary Regge 
trajectories, see details in |51 [T3]. 

3 Heavy baryon production within QGSM 

3.1 Sea charm and beauty quark distribution in proton 

Now let us analyze the charmed and beauty baryon production in the pp collision at LHC 
energies and very small p t within the soft QCD, e.g., the QGSM. This study can be interesting 
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for it may allow predictions for future LHC experiments like TOTEM and ATLAS and an 
opportunity to find new information on the distribution of sea charmed (c) and beauty (6) 
quarks at very low Q 2 . According to the QGSM, the distribution of c(c) quarks in the nth 
Pomeron chain (Fig. 1 (right)) is, see for example [T^] and references therein, 

f'$ ) ^)=C%5 c{ ^{l-x)^ (2) 

where a cn = —0^,(0), g cn — a p (0) — 2as(0) + (^(O) — (0)) + n — 1; (5 C ( E ) is the weight of 
charmed pairs in the quark sea, C^fs) ^ s ^ e normalization coefficient a^(0) is the intercept 
of the ip- Regge trajectory. Its value can be —2.18 assuming that this trajectory a^(t) is linear 
and the intercept and the slope a^(0) can be determined by drawing the trajectory through the 
J/\l/-meson mass mj/^ ~ 3.1 GeV and the x-meson mass m x = 3.554 GeV [TS]. Assuming that 
the ip- Regge trajectory is nonlinear one can get a^(0) ~ 0, which follows from perturbative 
QCD, as it was shown in 16J. The distribution of 6(6) quarks in the nth Pomeron chain 
(Fig.l(right)) has the similar form 

where a,b n = —atr(Q), gbn = °v(0) — 2«s(0) + (a p (0) — ar(0)) + n — 1; Q>(0) = 1/2 is the well 
known intercept of the p-trajectory; «b(0) ~ —0.5 is the intercept of the baryon trajectory, 
ax(0)) is the intercept of the T- Regge trajectory, its value also has an uncertainty. Assuming 
its linearity one can get a-r(O)) = —8, —16, while for nonlinear (bb) Regge trajectory a-r(O) — 0, 
see details in [T7| . Inserting these values to the form for f^(x) and f^(x) we get the large 
sensitivity for the c and b sea quark distributions in the nth Pomeron chain. Note that the 
FFs also depend on the parameters of these Regge trajectories. Therefore, the knowledge of the 
intercepts and slopes of the heavy-meson Regge trajectories is very important for the theoretical 
analysis of open charm and beauty production in hadron processes. 

Note that all the quark distributions obtained within the QGSM are different from the 
parton distributions obtained within the perturbative QCD which arc usually compared with 
the experimental data on the deep inelastic lepton scattering (DIS) off protons. To match these 
two kinds of quark distributions one can apply the procedure suggested in [TH] . The quantities 
9cn or gbn entering into Eq. ((J) and Eq. ([3]) arc replaced by the following new quantities depending 
on Q 2 

Q 2 + c" am am \~ • Q 2 + cl / 

The parameters c and d are chosen such that the structure function constructed from the 
valence and sea quark (antiquark) distributions in the proton should be the same as the one 
at the initial conditions at Q 2 = Qg for the perturbative QCD evolution. A similar procedure 
can be used to get the Q 2 dependence for the powers a cn and ctbn entering into Eqs. ()2|3j) [18]. 
Then using the DGLAP evolution equation [19] we obtain the structure functions at large Q 2 . 



gcn = ffcn(l + o2 ) ! 9bn = 9bn 1 + „ 2 ; ^ ) (4) 



3.2 Charmed and beauty baryon production in pp collision 

The information on the charmonium (cc) and botomonium (66) Regge trajectories can be found 
from the experimental data on the charmed and beauty baryon production in pp collisions at 
high energies. For example, Fig. 2 illustrates the sensitivity of the inclusive spectrum da/dx 
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of the produced charmed baryons A c to different values for a^(0). The solid line corresponds 
to a^(0) = 0, whereas the dashed curve corresponds to a^(0) = —2.18. Unfortunately the 
experimental data presented in Fig. 2 have big uncertainties; therefore, one cant extract the 
information on the a^(0) values from the existing experimental data. A high sensitivity of the 




Figure 2: The differential cross section da/dx for the inclusive process pp — > A C X at y/s 
62 GeV. 



inclusive spectrum da/dx of the produced beauty baryons Af, to different values for ax(0) is 
presented in Fig. 3 (left). The p t -inclusive spectrum of A;, has much lower sensitivity to this 
quantity, according to the results presented in Fig. 3 (right). Actually, our results presented in 
Fig. 3 could be considered as some predictions for future experiments at LHC, see Fig. 4. 

Now let us analyze the production of the beauty hyperon, namely A°, at small scattering 
angles 9^° ha the pp collision at LHC energies. This study would be reliable for the future 
forward experiments at LHC. The produced A" baryon can decay as A° — > J/^A , and J /if? 
decays into its branching ratio (Br = Tj/T) is 5.93 ± 0.06 percent, or into e + e~ (Br = 

5.93 ± 0.06%), whereas A can decay into pir- [Br = Tj/T = 63.9 ± 05%), or into nir° (Br = 
35.8 ±0.5%), see Fig.4. 

da f da 



d 3 pidM 3i J dM 12 dM 34 
<$ (3) (Pi +P2 -pi 2 )rfMi2 , 



(5) 



where 



da f 2 da pp ^A b x 

= / d p t A„ 



dM 12 dM u J 1 " dxd 2 p t A 

^3 



X 1/2 (Ml ff , M 2 2 , M 3 2 4 )A 1 /2(M 1 2 2 , M 2 , M^X^iM^, Af 2 , M 2 ) , 
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Figure 3: The differential cross section da/dx (left) and da/dP 2 (right) for the inclusive process 
pp -> A b X at V« = 4 TeV. 

Br^-j/* = (4.7 ±2.8) • 10~ 4 ; Brj/y^^ = (5.93 ± 0.06)%; Sr A o^ = (63.9 ±0.5)%. 
Here X(x 2 ,y 2 , z 2 ) = {{x 2 - (y + z) 2 )((x 2 - (y - z) 2 ) 
One can get the following relation 

d 3 Pi = ^pZpdfad^pdtp , (6) 

where £ p = Ap/p is the energy loss, t p = (j>i n — pi) 2 is the four- momentum transfer, 4>i is the 
azimuthal angle of the final proton with the three-momentum pi. 
Experimentally one can measure the differential cross section 

da 1 f da 

dtpdtpdMj,* = 2 P ^ P J d?p l dM M 4>1 { ' 

This distribution could be reliable for the TOTEM experiment, where J/^f decays into 

and A° decays into ir~p or for the ATLAS forward experiment, where A° decays as Ajj — > 

J/* A -> e+e~ ir°n (Fig.4). 

In Fig. 5 the distributions over Qj/^ (left) and £ p (right) are presented at different values 
of the intercept a-r(O) = (solid line), ax(0) = —8 (dashed line) and «x(0) = — 16 (dotted 
line), where 9j/^ is the scattering angle for the final J/ty. Fig. 5 shows a sensitivity of these 
distributions to the intercept of the ay Regge trajectory. Actually, the result presented in Fig. 5 
is a prediction for future LHC experiments on the heavy flavour baryon production at the LHC 
energies. 

4 Conclusion 

It was shown [6l [7] that the modified QGSM including the intrinsic longitudinal and transverse 
motion of quarks (antiquarks) and diquarks in colliding protons allowed us to describe rather 
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Figure 4: The decay A b — > J* A 1 



[i (e + e ) pir (nir ). 



satisfactorily the existing experimental data on inclusive spectra of heavy hadrons produced in 
pp and pp collisions It allows us to make some predictions for future LHC forward experiments on 
the beauty baryon production in pp collisions which can give us new information on the beauty 
quark distribution in the proton and very interesting information on the Regge trajectories of 
(bb) mesons. 
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